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Abstract. Data taken with the HEGRA system of imag- 
ing atmospheric Cherenkov telescopes during the 1997 
flares of Markarian 501 (Mkn 501) are reanalyzed using an 
algorithm providing improved energy resolution. A resolu- 
tion of 10% to 12% is obtained by accounting for the vari- 
ation of the Cherenkov light yield with the height of the 
shower maximum in the atmosphere. The improved en- 
ergy resolution is particularly relevant for the study of the 
high-energy cutoff in the spectrum, which might be caused 
by interactions with the intergalactic infrared background 
radiation. The reanalysis presented here confirms the re- 
sults obtained in the previous analysis, but hints a steeper 
slope of the spectrum in the region around 20 TeV. 

Key words: galaxies: BL Lacertae objects: individual: 
Mkn 501 - gamma rays: observations 



1. Introduction 

The attenuation of extragalactic gamma rays in the TeV 
energy range due to interactions with the cosmic infrared 
background radiation (Nikishov 1962; Gould & Schreder 
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1967; Stecker et al. 1992) is under intense discussion (see, 
e.g., Funk et al. 1998; Stecker & de Jaeger 1998; Aharonian 
et al. 1999a; Coppi and Aharonian 1999; Konopelko et al. 
1999a; Stecker 1999; Krawczynski et al. 2000; VassiHev 
2000; and references given therein), in particular since the 
energy spectra determined for the AGN Mkn 501 (Aharo- 
nian et al. 1999a, b,c; Samuelson et al. 1998; Djannati-Atai 
et al. 1999), obtained during its 1997 flares, show indica- 
tions of absorption features, particularly evident in the 
high-energy HEGRA data of Aharonian ct al. (1999a), in 
the following referred to as Paper I. The level of the cosmic 
infrared background is related to the evolution of galaxies 
and can serve to distinguish different cosmological models 
(see, e.g., Primack et al. 1999). However, due to the in- 
tense foreground radiation, a reliable direct measurement 
of the infrared background is a difficult task. Gamma rays 
in the TeV energy range provide an independent and com- 
plementary approach to study the infrared background 
radiation. Since the 77 — > e^e~ interaction cross section 
peaks near threshold, the absorption of gamma rays of en- 
ergy E probes the background level around a wavelength 
X/fim « nE/TeY, where k varies between 1 and 2 de- 
pending on the spectral shape of the background radia- 
tion. There is almost a one-to-one correspondence between 
the energy-dependent attenuation and the spectral density 
of the infrared background radiation. Recent analyses by 
Finkbeiner et al. (2000) of the background level at long 
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Fig. 1. (a) Attenuation of TeV gamma rays from Mkn 501 due to interaction with the infrared background radiation, 
for different assumptions concerning the density of the infrared background. Dotted hne: foUowing recent data from 
Finkbeiner et al. (2000), as shown by the sohd line of Fig. 1 in Protheroe & Meyer (2000). Dashed-dotted line: using 
the minimal flux consistent with the la statistical and systematic error bars of Finkbeiner et al. (2000). Full line: 
earlier LCFDM model of Primack et al. (1999). All curves use are calculated for the redshift z « 0.03 of Mkn 501, 
assuming a Hubble constant of 65 km s~^Mpc~^. (b) Spectral energy distribution E^dN/dE of TeV gamma rays from 
Mkn 501 at the source, reconstructed by dividing the measured distribution (HEGRA, Paper I, shown as triangles) 
by the attenuation. Filled circles: using the attenuation indicated by the dotted curve in part (a) of this Figure. Open 
circles: using the attenuation indicated by the dashed-dotted curve, shifting the HEGRA data down in energy by 
15% - the maximal systematic error on the energy scale - and including the systematic error bars of the HEGRA 
measurement. 



wavelengths, at 60 fiia and 100 /im, show a significantly 
higher flux than expected, which results in strong absorp- 
tion of Mkn 501 gamma rays above 10 TeV (Fig. |(a)). 
While earlier analyses using a lower level of the infrared 
background (e.g., the models of Primack et al. (1999), also 
illustrated in Fig. |l|(a)) found good consistency with the 
TeV data (see, e.g., Konopelko et al. (1999a)), these new 
results appear to be in contradiction with measurements 
using the HEGRA Cherenkov telescopes, which show a 
gamma-ray flux from Mkn 501 up to at least 16 TeV 
(Paper 1). In order to explain the observed flux in the 
presence of strong high-energy attenuation, the spectral 
distribution of gamma-rays at the source would have to 
exhibit a rise at high energies (Fig. 0(b)), inconsistent 
both with acceleration models and with energetics of the 
source (see, e.g.. Paper I and Coppi and Aharonian 1999; 
Guy et al. 2000; Protheroe & Meyer 2000; Finkbeiner et 
al. 2000; and references given therein). Exotic solutions 
of this "gamma-ray crisis" have been proposed, such as 
violations of Lorentz invariance (Kifune 1999) or photon 
condensates (Harwit et al. 1999); the latter explanation 
has been questioned on rather general grounds (Levin- 
son 2000) and seems already to be ruled out by a recent 
HEGRA analysis (Aharonian et al. 2000). 



The TeV gamma-ray data and the infrared background 
data can be made marginally consistent by fully exploiting 
the systematic errors of both measurements, i.e., taking 
the minimal infrared flux values and gamma ray flux val- 
ues and reducing the HEGRA energy scale by 15%. In this 
case, a marginally acceptable source spectrum is achieved, 
see Fig. |^(b), in particular given that the highest HEGRA 
data point is, at the 2a level, consistent with spill-over 
from the lower energy bins. The bulk of the change in the 
calculated source spectrum - a factor 100 for the highest 
data point - can be attributed to the uncertainty in the 
infrared data; the effect of the HEGRA uncertainties is 
smaller - about a factor 10 for the highest point. 

The interpretation of the experimental results is cur- 
rently limited by the systematic errors both of the mea- 
surement of the infrared background, and of the TeV spec- 
tra. One potentially very important factor is the energy 
resolution in the detection of the TeV gamma rays. The 
HEGRA system of Cherenkov telescopes, which provides 
the highest-energy data, has an energy resolution around 
20% (see Aharonian et al. 1999b, Konopelko et al. 1999b, 
Paper I) . While this resolution is ample for the reconstruc- 
tion of the usual power-law energy spectra, it may smear 
out the extremely sharp cutoff features expected for ab- 
sorption due to infrared background, in particular for high 
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long- wavelength infrared background levels (Fig. 0(a)). 
The highest-energy gamma rays detected could simply 
represent spill-over from lower energies due to the finite 
energy resolution. Indeed, while Aharonian et al. (Paper 
I) measure photons with reconstructed energies up to 21.5 
TeV and beyond, they state that a sharp cutoff at 16 TeV 
cannot be excluded by the data. It is obviously of con- 
siderable interest to improve the energy reconstruction of 
the system of telescopes, in order to more precisely map a 
possible cutoff. A scheme to improve the energy resolution 
of systems of Cherenkov telescopes by roughly a factor 2 
has been discussed by Hofmann et al. (2000), making use 
of the abundant information provided by such systems 
for each individual gamma-ray air shower. In this note, 
we present a reanalysis of the HEGRA Mkn 501 spectral 
data with such improved analysis algorithms. 

2. Determination of gamma ray spectra 

One of the key limiting factors in the energy determina- 
tion of air showers with imaging Cherenkov telescopes is 
the fluctuation in the height of the shower maximum in 
the atmosphere. Deeper showers have a larger light yield 
both because of the reduced distance to the telescopes and 
because of the reduction of the threshold for Cherenkov 
light emission with decreasing height, which allows more 
particles to contribute to the emission of Cherenkov light. 
As demonstrated by Hofmann et al. (2000) and by Aharo- 
nian et al. (2000), the multiple views of a shower provided 
by the HEGRA stereoscopic system of Cherenkov tele- 
scopes allow to determine the height of the shower max- 
imum with an error of 500 m to 600 m. Using look-up 
tables which depend both on the distance to the shower 
core and the shower height to relate the Cherenkov light 
yield to the gamma-ray energy, the energy resolution of 
the HEGRA system can be improved significantly (Hof- 
mann et al. 2000). A further improvement can be reached, 
if the known location of the gamma-ray source is used as 
input in the reconstruction of the shower geometry, pro- 
viding a greatly improved localization of the shower core, 
and hence a better measurement of the distance of a given 
telescope to the shower core. Combining all these tech- 
niques, an energy resolution around 10% to 12% is ob- 
tained (Hofmann et al. 2000), implying - for the identical 
set of events - a factor 2 improvement over more conven- 
tional techniques such as used in Aharonian et al. (1999b) 
and in Paper I. 

The new analysis technique was applied to the data 
collected with the HEGRA telescope system during the 
extraordinary 1997 flares of Mkn 501. The telescope hard- 
ware, the event samples, the event reconstruction and 
event selection, and the determination of energy spectra 
are described in Aharonian et al. (1999b), and in Paper I. 
In particular, relatively loose cuts on the shower direction 
and in the image shapes are used to select gamma rays 
from Mkn 501, in order to minimize the efficiency correc- 



tions and the systematic errors resulting from these cor- 
rections, and events are only used up to core distances of 
200 m from the central telescope. In the current analysis, 
the new procedure is employed to determine the shower 
energy on the basis of the measured height of the shower 
maximum, and the core location is determined given the 
known source location. We will, in the following, only ad- 
dress the energy range above 3 TeV, where the effective 
detection area of the telescope system is approaching a 
constant and where threshold effects are negligible. The 
look-up tables used in the energy determination were op- 
timized for this energy range. At lower energies, errors in 
the determination of the energy spectra are dominated by 
the modeling of trigger efficiencies (see Paper I), and no 
attempt was made to refine the spectra in this range. Ac- 
cording to Monte Carlo simulations, the improved energy 
estimate provides an energy resolution of 10% to 12% and 
has a systematic deviation between the true energy and 
the average reconstructed energy of less than ±3%, in the 
range between 3 TeV and 30 TeV. Fig. || compares the 
performance of this energy reconstruction with the pre- 
vious technique; the new approach provides both better 
resolution and reduced non-Gaussian tails. 
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Fig. 2. Ratio of the reconstructed and true gamma-ray 
energies, as derived from Monte-Carlo simulations at 20° 
zenith angle, for energies between 3 TeV and 30 TeV. His- 
togram: new algorithm, points: conventional analysis. 



3. Tests of the reconstruction algorithm 

A key issue is of course whether the Monte-Carlo simu- 
lation (Konopelko et al. 1999b) - from which the energy 
look-up tables are derived - is accurate enough to guaran- 
tee that the identical energy resolution is obtained for the 
simulations and the real data. Lacking a monoenergetic 
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test beam, the energy resolution of Cherenkov telescopes 
cannot be determined directly. Fortunately, however, the 
information provided by the up to four views of a given 
shower is highly redundant, and allows a number of cross 
checks, some of which are briefly summarized in the fol- 
lowing. We address in sequence the reconstruction of the 
shower core, the determination of the shower height, and 
the energy determination and influence of the atmosphere. 

In most of the events in the range above 3 TeV, all four 
telescopes trigger and provide images; the mean number of 
telescopes used in the reconstruction is 3.7. The precision 
in the location of the shower core can be tested by recon- 
structing, in four-telescope events, the same shower inde- 
pendently by two stereo subsystems of two telescopes each. 
The rms difference in the x or y coordinates of the core, 
comparing two telescopes with the other two, is about 
3.1 m for small core distances, and grows to 5.6 m for dis- 
tances of 100 m. These differences are about 30% larger 
than predicted by the simulation, indicating remaining 
alignment errors. The same simulation predicts an error 
in the core coordinates of less than 2 m rms over the full 
distance range, when combining all telescopes. Scaled up 
by 30%, this error is, however, still safely below values 
which could affect the energy determination. 

For the energy determination, it is crucial that the 
simulation provides a proper modeling of the atmosphere, 
and in particular that the average height of the shower 
maximum is reproduced. Otherwise, the energy correc- 
tion based on the height of the shower maximum might 
introduce a serious bias. Fig. ^ shows, as a function of 
zenith angle, the measured shower height, compared with 
simulations. In order not to introduce a bias by cuts on 
the - potentially biased - shower energy, events were se- 
lected by requiring two telescopes with at least 100 pho- 
toelectrons. Data and simulation agree to better than 2%. 
Simulations show that even a 5% mismatch in the aver- 
age shower height has an almost negligible effect on the 
energy resolution, and introduces a bias in the overall en- 
ergy scale of less than 3-4%. Also the widths of the dis- 
tributions of shower heights agree within statistical errors 
between data and simulation; at 20° zenith angle, e.g., the 
measured (Gaussian) width is 0.82 ± 0.02 km, compared 
to 0.80 ± 0.01 km in the simulation. We note that due 
to the selection condition used for the data of Fig. |[ the 
energy of gamma rays increases with increasing zenith an- 
gle. The zenith-angle dependence of the average shower 
height is caused by a combination of two effects: at larger 
angles (and fixed energy), the shower maximum moves up 
since the amount of radiation lengths traversed at a given 
height increases like 1/ cos(^2enith)' however, showers at 
larger zenith angles require a higher energy to provide two 
images above 100 photoelectrons and therefore penetrate 
deeper. 

Checks of the energy resolution use the feature that 
for a given shower and a core location determined on the 
basis of the orientation of the images, the light yield ob- 
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served in each individual telescope can be converted into 
an energy estimate. For the final energy value these esti- 
mates are averaged, but the consistency of the up to four 
energy values for each shower provides stringent tests of 
the quality of the reconstruction procedure. 

A crucial input in the energy reconstruction is the de- 
pendence of the measured light yield on the distance to the 
shower core. This input can be checked by comparing the 
energy estimates of two telescopes located at different dis- 
tances from the shower core. We find that the (averaged) 
ratio of the energy estimates provided by two telescopes 
shows a systematic variation with the two core distances 
(in the range between and 200 m), at a level of 4.2% rms. 
A 4.2% systematic variation should not spoil the 10% en- 
ergy resolution; furthermore, the effect is fully reproduced 
by the simulations, which show, under the same condi- 
tions, a 4.4% rms variation. 

The previous test checks if the simulations reproduce 
the variation of the mean light yield with core distance. 
An equally important test is to see if the fluctuations are 
properly described. Using four-telescope events, the me- 
dian rms spread of the four individual energy estimates 
was determined to 12.3 ± 0.3%, slightly below the Monte 
Carlo value of 13.1%. The conventional energy reconstruc- 
tion shows a 19.5 ± 0.5% median rms spread between tele- 
scopes, compared to the simulation value of 19.2%; shower 
height fluctuations result in strong correlations between 
the individual telescopes, and govern the energy resolu- 
tion. 
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Finally, one has to worry about varying atmospheric 
conditions, which might introduce time-dependent varia- 
tions in the light yield. Since such effects will most likely 
influence all telescopes in more or less the identical man- 
ner, they would not show up in inter-telescope compar- 
isons. One can distinguish two effects: long-term seasonal 
variations and short-term variations. The influence of sea- 
sonal variations on the Cherenkov light yield is discussed 
by Bernlohr (2000). Atmospheric density profiles and tem- 
perature profiles on the Canary Islands are continuously 
monitored by radiosondes; on the basis of these data 
(Cuevas 1997) and of the simulations of Bernlohr (2000) 
we conclude that the seasonal variation of the Cherenkov 
light yield should be less than 2% rms for the Mkn 501 
data set. More critical are short-term variations in atmo- 
spheric transmission. As discussed in Paper I, the data set 
used here is selected to exclude nights with excessive at- 
tenuation. To monitor the attenuation, the trigger rate of 
the telescope system is used; the rate varies roughly like 
the 1.6th power of the energy threshold and variations in 
the transmission cause corresponding rate variations, am- 
plified by a factor 1.6. For the data set after selection, the 
mean trigger rate (at a given zenith angle) varies day-by- 
day by 4% rms for the first period of data, and by 6% 
rms for periods 2 and 3 (Aharonian et al. 1999b; by mis- 
take, the paper quotes 4% FWHM instead of the correct 
value of 4% rms). These rate fluctuations translate into 
fluctuations in the light yield - and hence of the energy 
scale - of 3% to 4% rms, sufficiently small compared to 
the resolution of 10% to 12% rms. 

On the basis of these and other tests, we believe that 
in the data there are no systematic effects or additional 
fluctuations, which would significantly worsen the energy 
resolution compared to the resolution predicted by the 
simulations. 

Concerning the absolute energy scale, the new anal- 
ysis is exactly like the previous analyses susceptible to 
imperfections in the modeling of the atmosphere, of the 
transmission of the telescope optical systems, and of the 
photomultiplier and electronics gains. Therefore, the 15% 
uncertainty on the energy scale discussed in Paper I also 
applies here. The same is true for systematic errors at high 
energies due to photomultiplier nonlinearities or electron- 
ics saturation. 

4. Results and interpretation 

Fig. 1^ shows the spectral energy distribution obtained 
with the new energy reconstruction, together with the pre- 
vious results of Paper I. The new data points are listed 
in Tab. |^. Except for the highest point in the previous 
data set, at E = 21.45 TeV, the two analyses are in very 
good agreement. In Paper I, the influence of the energy 
resolution was addressed by studying a model, where the 
source spectrum shows a power law with an index around 
2, an exponential cutoff, and in addition a sharp cutoff at 
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Spectral energy distribution E'^dN/dE in units 
erg cm~^s~^. See text for systematic errors. 



an energy £^cut- the 2a level, a hard cutoff below 16 
TeV was excluded in Paper I; this result implies that at 
the 2(7 level, in particular the data point at 21.45 TeV is 
consistent with spill-over from lower energies. In the new 
analysis, no significant signal is observed in this energy 
range, and upper limits are below the value of the previ- 
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ous point. We conclude therefore that at these very highest 
energies, the spectrum appears indeed steeper than indi- 
cated by the data of Paper I, but consistent within the 
systematics quoted in Paper I. While this change goes in 
the expected direction, a statistical fluctuation at the 2-3ct 
level cannot be excluded as the origin of the change, due 
to the low statistics for this data range in either analysis. 
In the range up to the data point at 17.08 TeV of Paper 
I, there is almost no change observed due to the improved 
energy resolution. Combining the previous data below 3 
TeV with the new data points above 3 TeV, a good fit to 
the data is achieved with the parametrization 



dN/dE - E 



-2.03±0.03 



e^p(-{E/Eo 



1.31±0.11 



with Eo = 8.21 ± 0.63 TeV. The fit indicates that indeed 
a super-exponential cutoff is favored. We note - as al- 
ready mentioned in Paper I - that the fit parameters are 
strongly correlated and cannot be varied independently 
within their errors. Errors quoted for the fit parameters 
are purely statistical. 

Carrying out a similar analysis with a hard cutoff for 
the new results, a cutoff below 16.5 TeV can be excluded. 

In summary, the conclusions concerning the level of 
the infrared background allowed by the HEGRA Mkn 501 
energy spectra remain unchanged: taken at face value, the 
HEGRA data and the infrared flux levels of Finkbeiner et 
al. (2000) are inconsistent; making full use of the system- 
atic errors allowed for the HEGRA data and in particular 
for the infrared flux determinations, the two data sets can 
be brought to a marginal consistency, without a clear need 
to invoke exotic explanations. 
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